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University of British Columbia 

We consider a multitype branching process with immigration in 
a random environment introduced by Key in [Ann. Probab. 15 (1987) 
344-353]. It was shown by Key that, under the assumptions made in 
[Ann. Probab. 15 (1987) 344-353], the branching process is subcritical 
in the sense that it converges to a proper limit law. We complement 
this result by a strong law of large numbers and a central limit the- 
orem for the partial sums of the process. In addition, we study the 
asymptotic behavior of oscillations of the branching process, that is, 
of the random segments between successive times when the extinc- 
tion occurs and the process starts again with the next wave of the 
immigration. 

1. Introduction and statement of results. In this paper we consider the 
multitype branching process with immigration in a random environment 
(MBPIRE) introduced by Key in [14]. Broadly speaking, the MBPIRE is 

a vector- valued random process Z n = {Z n l \ Z n 2 \ . . . , Zn) that describes 
evolution of a population of particles of d different types in discrete time 

(i) 

n = 0,1, The integer number d > 1 is fixed and Z n denotes the num- 
ber of particles of type i in generation n. We next present the process 
Z = (Z n ) n >Q in more detail. 

The environment lo = (w ri ) n gz is a stationary ergodic sequence of random 
variables uo n taking values in a measurable space (5, B). We denote by P the 
law of uj on (SI, ^*) := (S Z ,B® 1 ') and by Ep the expectation with respect to 
P. Both the law of the immigration and the branching mechanism depend 
on the realization of the environment. 

Let Z + denote NU {0}. The immigration X = (X n ) n€ % + is a sequence of 

independent conditionally on u random vectors X n = (X n l \ . . . , ) £ 
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describing the number of immigrants of all types arrived at time n G Z + . 
Let Vd be the set of probability measures on Z^ and assume that the law 
Qu;,n of X n in a fixed environment to is a function of the past state u> n ~i of 
the environment: 

(1-1) = g(^n-l) 

for some measurable function q:S -^Vd- Here we regard Vd as a subspace 
of the space of bounded real-valued sequences endowed with the Borel 
cr-algebra induced by the product topology. 

Every particle in the branching process is a descendant of an immigrant 
but the immigrants themselves are not counted in the population. At each 
unit of time, every particle of type i present in the system (including the 
immigrants) splits (for a given realization of the environment independently 
of the previous history and of the other particles) into a random number of 
offspring of all types according to a law p$ jTl which is a function of uj n : 

(1-2) P^n-PiM 

for some measurable functions pi\S ^Vd,i = 1, • • • , d. 

Thus, conditionally on the environment, (Z n ) n ^i + is a nonhomogeneous 
Markov chain that satisfies the initial condition Zq = [here and throughout 
stands for (0, . . . , 0) G Zl] and the branching equation 

(1.3) Z n+1 = J2 E L Sm forn>0, 

i=l m=l 

where the random vectors Ln]m = (Lhjm £n,m ) G Z^ (with n G Z + ; i = 
1, . . .,d;m G N) are independent and Ln,m are distributed according to the 
law pi^n for every m G N. 

The underlying probability space (fi x T,^"® £,P), where T is a space of 
family trees that describes the "genealogy" of the particles and Q is its Borel 
cr-algebra, is constructed using the recipe given in [12], Chapter VI. The mea- 
sure P on x T is defined as P <8> -FL>, where for co G f2, -FL is the (popularly 
known as quenched) law on (T,Q) consistent with the preceding description 
of the process in the environment uj. The marginal P(-) = Ep(P aJ (-)) of P 
on (T,Q) is referred to as the annealed law of the process. The expecta- 
tions with respect to P^ and P are denoted by E w and E, respectively. For a 
random vector v = (v^\ . . . ,v^) G M d , Ep(v) [correspondingly E u (y), E(i>)] 
stands for the vector whose ith component is Ep(v^)(E UJ (v^), E(«W)). 

To formulate our assumptions on the process (Z n ) n£ z + we need to in- 
troduce some additional notation. For any constant > and real-valued 
random variable X we set 

(1-4) \\X\\ u , p :=(E u (\Xf))W. 
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We often use L n l as "a random vector distributed according to the law 
Pu>,n" and in such cases omit the lower index 1 writing it simply as Ln ■ The 

(i) (i j) 

jth component of the vector L n is denoted by L n . For any n € Z + and 
i, j = 1, . . . , d we define the following functions of the environment u: 

1 n,/3\ i0 ) ■— W^n IL,/3> lvl n,f3 V^J ■— W^n IL,/3- 

We denote by J nj( g the vector (I n g, ■ ■ ■ , I n d ^) and (with a slight abuse of 
notation) by M n ^ the matrix whose ith column (not the row) is the vector 
(M^f,. . . , M^jf ). To simplify the notation we write 

In ■= I n ,i and M n := M n>1 . 

Finally, for v = (v^, . . .,v&) G K d we set \\v\\p =: (Ef=i |w W | /3 ) 1//3 , and for 
a d x d matrix A denote by \\A\\p the corresponding operator norm. 
The following is the basic set of conditions imposed in this paper. 



Assumption 1.1. We have: 

(Al) .Ep(log + ||/o||i) < °°> where log + x := max{0, logx}. 

(A2) £7 P (log + ||M ||i)<oo. 

(A3) lim n ^ 00 n- 1 £;p(log||M n _i---MiMo||i)<0. 



Note that (Al) implies P(ljp < oo) = 1 and P(M^' j) < oo) = 1 for all 
and that (A2) yields the existence of the limit in (A3) by the sub-additive 
ergodic theorem of [11]. 

The model of MBPIRE was introduced by Key in [14], where it is also 
shown that the process is sub critical under Assumption 1.1, that is, Z n 
converges in distribution to a proper limit law (see Theorem 2.1 below where 
a more precise statement is cited from [14]). 

The model is a generalization on one side for the multitype branching pro- 
cesses in random environment (without immigration) introduced by Athreya 
and Karlin in [3] and on the other side for a special single-type case con- 
sidered by Kesten, Kozlov and Spitzer in [13]. The latter process is closely 
related to one-dimensional random walks in random environment (RWRE) 
and from this point of view was studied by many authors (see, e.g., [1, 8, 13] 
and Section 2 of the survey [17]). The limiting distribution of Z n [tt w de- 
fined below in (4.1) turns out to be a geometric law in this case], the weak 
law of large numbers and the central limit theorem for the partial sums of 
this process in a stationary ergodic environment can be obtained by using 
corresponding statements for the associated RWRE. 

In this paper we study the asymptotic behavior of the sequence (Z n )nez + 
exploiting both tools developed in the theory of random motion in random 
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media (e.g., the method of the "environment seen from the position of the 
particle") as well as the construction of the limiting distribution for the 
general MBPIRE provided in [14]. 

The general idea of the proofs is as follows. Let Ui n G M. d represent the 
progeny in generation n of all immigrants arrived at time i. Then Z n = 
J2l=o Ui,n stochastically grows to its stationary distribution as n goes to in- 
finity. The limiting law corresponds to the random variable Z n = J27=-oo U% n 
which is a.s. finite under Assumption 1.1. Here we introduce "demo" immi- 
grants arriving at negative times and keep the notation Ui n for their progeny. 
The process Z n is stationary and moreover, it turns out that it inherits some 
ergodic and mixing properties of the environment. Furthermore, since the 
MPBREs (Ut 

n)n>i extinct a.s. for any i G Z (cf. [3]), the random variables 
Z n and Z n coincide for all n large enough. Therefore one can obtain limit 
theorems for the process Z n from standard results for the ergodic (mixing 
under additional assumptions) stationary process Z n . Note that in this as- 
pect the process with immigration is quite different from the supercritical 
MPBRE without immigration, where Z n grows in exponential rate (cf. [6]). 
We next turn to the presentation of our main results. Let 

n-1 

Sri = y ] Z n . 

The following strong law of large numbers is proved in Section 2: 

Theorem 1.2. Let Assumption 1.1 hold. Then P(lim n ^ 00 ^ = p) = 1, 
where 

oo 

(1.5) p := J2 E P (M n ■ ■ ■ MiMq/o) € (K+ U {cx)}) d . 

n=0 

The proof of the theorem is by coupling (Z n ) n£ z + with a stationary and 
ergodic process (Z n ) n( zz [introduced below in (2.2)] such that Z n are dis- 
tributed according to the limit law of Z n , and applying then the ergodic 
theorem to (Z n ) ne z (cf. Lemma 2.2). 

In the case where (w n )nez is a strongly mixing sequence with exponential 
rate we obtain in Section 3 the following central limit theorem. Let 

(1.6) a(n) = sup{\P(A n B) - P(A)P(B)\ : A G T n \ B G ^ }, 
where, for n G Z, 

(1.7) J^ 1 := aiio-i :i>n) and T n := aiuJi : i < n). 

Recall that the sequence ui n is called strongly mixing if ol{ti) — > n »oo 0. For 

strongly mixing sequence it holds that (cf. [9], page 10) 

(1.8) \E P (fg) - E P (f)E P (g)\ < 4a(n) 



BRANCHING PROCESSES IN A RANDOM ENVIRONMENT 



5 



for all ^"-measurable random variables f{oo) and all J-Q-va.easxaab\e random 
variables g{ui) such that P(\f \ < 1 and \g\ < 1) = 1. 
For any n G Z + and i G {1, . . . , d} let 

:= ^{progeny of type i over all generations 

of all X n immigrants arrived at time n} 

and denote Y n := (Y^\ . . . , Y^ ) G Z^. An equivalent definition of Y n is 
given in (3.1) below. 

Theorem 1.3. Let Assumption 1.1 ZioW and suppose in addition that: 

(Bl) limsup n ^ 00 n^ 1 log£;p(||M n _i iK ---Mo iK /o, K ||K) <° f or some K > 2 - 
(B2) limsup n ^ OC) n~ l loga(n) < 0, where the mixing coefficients a(n) are 
defined in (1.6). 

Then n _1 / 2 (S' n — n- p) converges in distribution to a Gaussian random vector 
Sqo with zero mean. Moreover, if the following condition holds for some 
ie{l,...,d}, 

(B3^) There is no function fi-.n^Z^ such that P^Y^ = fi(u)) = 1 
for P-a.e. environment u, 

then EQSS] 2 ) > for that i. 

Remark 1.4. (i) Moment condition (Bl) is used in the proof of The- 
orem 1.3 to ensure that E(||Yo||k) < 00 ■ Mimicking the estimates on the 
moments of a single-type branching process carried out in the course of the 
proof of [8], Lemma 2.4 (see also [17], Lemma 2.4.16), it is not hard to check 
that the conclusion of Theorem 1.3 still holds if (Bl) is replaced by 

(Bl') There exist (nonrandom) constants k > 2 and m > such that 

P(||Xo||o;,k < m) = 1, P(||Lq \\ w<k < m) = 1 for all i = 1, . . . ,d, and in ad- 
dition limsup n ^ 00 n- 1 log J Bp(||M n _i---Mo||^) < 0. 

Note that (Bl') implies by Jensen inequality that the conditions of The- 
orem 1.2 hold. 

(ii) The random vectors Y n , n G Z+, are independent in a fixed envi- 
ronment, and one can check that if Assumption 1.1 holds and condition 
(Bl) is satisfied, then the Lindeberg condition (cf. [10], page 116) is fulfilled 
for the sequence (Y n ■ i) n ez + for every t ^ G M. d . Therefore, for P-a.e. u>, 
the sequence (Y n ■ i) n ez + obeys a CLT with a random centering under the 
quenched measure P w . Alili (cf. [1], Theorem 5.1 and Section 6) provided a 
set of conditions on the environment [excluding i.i.d. sequences (w n )nez + ] 
sufficient to replace the random centering in the quenched CLT of this kind 
by a constant. 
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This theorem is similar to the CLT for hitting times of one-dimensional 
random walks in a stationary ergodic environment (cf. [17], see also [2], 
Theorem 4.3) that can be equivalently restated in terms of the single-type 
process considered in [13]. The analogy between hitting times of RWRE and 
random variables Y n that exists in the special case can be carried over, to 
some extent, to the general MBPIRE. Again, the limiting distribution of Z n 
is an important ingredient of the proof. 

Our next result deals with the regeneration times u n defined by 

(1.9) ^o = and v n = inf{i > v n _\ : Zi = 0}, 

with the usual convention that the infimum over an empty set is oo. 

Theorem 1.5. Let Assumption 1.1 hold and suppose in addition that 
there exists constant m G N such that 

(1.10) P(P w (Z m = 0|Z = 1)>0)>0, 

where 1 := (1, . . . , 1) G R d . 

Then F(Z n = i.o.) = 1. Moreover, if 

(1.11) P(P U (Z 1 =0\Z = 1)>0) = 1, 
then P(lim n _ >00 v n /n = fx) = 1 for some constant fx > 0. 

The proof of this theorem given in Section 4.1 uses the same coupling with 
the sequence Z n as in the proof of Theorem 1.2 along with a change of mea- 
sure argument. Assumption (1.10) is essentially a condition of Key implying 
that the limiting distribution of Z n puts a positive probability on (see 
Lemma 4.2 below). 

The law of large numbers for v n is derived from a somewhat more gen- 
eral result stated in Proposition 4.3. Let Zua denote the segment of the 
branching process between times i and j, that is, 

(1.12) Z[ i:j y.= (Zi,Z i+1 ,...,Zj), i<j, 

and let 6 denote the shift operator on the space of environments f2, that is, 

(1.13) (6u) n = u n+ i. 

We introduce in Section 4 a probability measure P on the underlying mea- 
surable space (ft x T, J 7 x Q) such that under P the sequence of triples 
{{0 Vn - 1 uj,v n - i/ n -i, Z[ Un _ 1+liUn ])) ne n is stationary and ergodic. The distri- 
bution of the sequence (Z n ) n >Q under P is the Palm measure P((Z n ) n >o G 
•) = F((Z n ) n >Q £ -\Zq = 0) and it turns out to be equivalent to P((Z n ) n >o G •) 
if (1.11) holds. 
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We next consider the case where the law of the branching process is 
"uniformly elliptic" and the environment uj is a uniformly mixing sequence. 
In this case we show that the sequence of pairs {{v n — u n -\, Z^ 1/n _ 1+ i L ,])) n ^ 
is uniformly mixing and converges in law to its stationary distribution. Let 

<p(n) = sup{\P(A\B) - P(A)\ lief.Be Fq,P{B) > 0}, 

where the cr-algebras T n and J-q are defined in (1.7). Recall that the sta- 
tionary sequence (w n ) ng z is called uniformly mixing if (p{n) — ^n— >oo 0. For 
uniformly mixing sequence it holds that (cf. [9], page 9) 

(1.14) \E P (fg) - Ep(f)E P (g)\ < 2<p(n)E P (\f\) 

for all ^"-measurable random variables f(to) and all .Fo-measurable random 
variables g(u) such that P(\f \ < 1 and \g\ < 1) = 1 [cf. with (1.8)]. 
The following theorem is proved in Section 4.2. 

Theorem 1.6. Let Assumption 1.1 hold and suppose in addition that 
the sequence {uS) n ^% is uniformly mixing and that there exist d+1 measures 
e W £ ~p dj z = l,...,<i + 1, and constant e S (0, 1) such that the following 
hold: 

(i) For every v G lA, 

p (P%( v ) > ^(v) fori = l,...,d and q^ (v) > ee^ d+1 \v)) = 1, 

(ii) eW(0) >0 for all i = 1, . . . ,d. 
Then: 

(a) The sequence {(v n — v n -i, ^[ 1 /„_ 1 +i, i /„]))neN is uniformly mixing. More 
precisely, 

sup sup sup dP^I^-P^D-^n^O, 

meN Asg^+m B£g„ m ,p(b)>o 

where we denote Q n = a(Z{ :i>n) and Q n = a(Zi,uJi :i <n). 

(b) The sequence ({v n — v n -\, ■Z'[^_ 1 +i,i/ ?l ]))neN converges to its stationary 
distribution. 

Note that the assumptions of the theorem imply that condition (1.11) of 
Theorem 1.5 is valid. 

The rest of the paper is organized as follows. Section 2 contains the proof 
of the law of large numbers (Theorem 1.2), Section 3 is devoted to the proof 
of the central limit theorem (Theorem 1.3), and the proofs of Theorems 1.5 
and 1.6 are included in Section 4. 
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2. The LLN: Proof of Theorem 1.2. Theorem 1.2 is obtained in this 
section from the corresponding statement for the stationary and ergodic 
"version" (Z) neZ of the process (^n)n£Z_|- which is introduced below in (2.2). 

For any i G Z + and n > i define the vector Ui jn = (U^, . . . , U^) by set- 
ting Ui t i = Xi and, for n> i, 

= ^{progeny of type j in generation n 

(2.1) 

of all immigrants arrived at time i}. 

Thus Z n = J21=o Ui, n for n G Z+. The sequence Hi := (Z7j jn ) n >j form a mul- 
titype branching process in random environment as introduced by Athreya 
and Karlin in [3] (abbreviated as MBPRE in what follows). Assume that 
the underlying probability space is enlarged to include random vectors Xi 
and processes Hi for i < such that: 

(i) For every i G Z, Uu = Xi. 

(ii) The processes {Ui)i^% are independent in a fixed environment. 

(iii) For any i < 0, PujiUi G •) = Pgi w (Uo G •), where the shift is defined 
in (1.13). 

That is, (Ui)i£i is a stationary collection of MBPREs with independent lines 
of descent in a fixed environment, each Hi is starting from the immigration 
wave described by Xj. 
Let 

71-1 

(2.2) Z n = ^ U i>n , nGZ. 

i=— oo 

Some components of the vectors Z n may be infinite a priori. However, the 
following was in fact proved by Key (cf. [14], Theorem 3.3; the random vari- 
ables Z n are not defined explicitly in [14] but they are "present implicitly," 
e.g., in Lemma 2.2 there). 

Theorem 2.1 ([14]). Let Assumption 1.1 hold. Then Zq is a proper 

random vector, that is, W(Zq < oo Vz = 1, . . . , d) = 1, and ¥(Zq = v) = 
linin^oo P(Z n = v) for every v G Z^_. 

That is, (Z n ) ng z provides a "stationary version" of the process Z n , and 
the latter converges in law to its stationary distribution. We next show that 
the sequence (Z n ) n ^z is ergodic and derive the law of large numbers for the 
partial sums of Z n from this fact. 

Lemma 2.2. Let Assumption 1.1 hold. Then: 
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(a) The sequence U = (U n ) n€ z is ergodic. 

(b) (Z n ) n€ z is a stationary ergodic sequence. 

(c) P(3 n G N : Z n = Z n /or n > n ) = 1. 

(d) For any measurable function /:Z+ — ► M i/ie following holds P-a.s., 
provided that the expectation in the right-hand side exists: - ^27=0 fi^i) — >n-> 
E(/(Z )). 

Proof, (a) Let ^4 G a(U n :n G Z) be an invariant set of the shift op- 
erator defined by (0U) n =U n+ i, that is, A = -1 A modulo a P-null 
set. Then P-a.s., P U (A) = P U) {®~ 1 A) = P e -i w (A). Therefore the function 
f(uj) := P W {A) is invariant under the ergodic shift 6 on the sequence to = 
(u> n ) n( zz, that is, f{uj) = f(6co), P-a.s. Hence f{uj) is a P-a.s. constant func- 
tion. Since the random variables lA n are independent under P w , Kolmogorov's 
0-1 law implies that f(u) G {0, 1}. It follows that F(A) = E P (f(u)) G {0, 1}. 

(a) The claim follows from part (a) of the lemma and the definition (2.2) 
of Z n . 

(c) Follows from the extinction criterion for MBPRE given in [3], The- 
orem 12 (with notation of Remark following this theorem). The criterion 
is applied to the progeny of all Zq particles living at time in the process 
{Z n ) n& i, and implies that the extinction of the branching process formed by 
these particles occurs eventually under Assumption 1.1. 

(d) Follows from the ergodic theorem applied to the sequence (Z n ) n£ x 
and part (c) of the lemma. □ 

Part (d) of the lemma contains the claim of Theorem 1.2 and we now turn 
to the proof of the CLT for the partial sums of Z n . 

3. The CLT: Proof of Theorem 1.3. The CLT for the process Z n is de- 
rived here from the corresponding statement for the sequence Y n (introduced 
right before the statement of Theorem 1.6). The definition of Y n can be 
equivalently written as follows: 

oo 

(3.1) Y n = u n,i, neZ, 

i=n+l 

where the random vectors U n i are introduced in (2.1). It is not hard to 
check that each component of the vector E u (Yq) = Y^=i Mi-i " ' ' MqIq is P- 
a.s. finite under Assumption 1.1 (see, e.g., [14], Lemma 3.1). Therefore Yq is 
a proper random vector, E(Yo) = P, an d it follows from Lemma 2.2(a) that 
(Y n ) n< zz is a stationary ergodic sequence. In order to prove the CLT for the 
partial sums of Y n we shall use the following general CLT for mixing vector- 
valued sequences. Recall the definition of the norm \\v\\ K = (J2i=i \v^\ K ) l ' K , 
v£R d . 
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Lemma 3.1. Let (Y n ) n& z + be a stationary sequence of random vectors 

Y n = (T n (1) ,...,y n (d) ) in R d such that E(||F ||«) < oo for some k>2. For 
n £ N let T n = o~(Y m :m>n) and denote 

(3.2) Xn = sup{|P(^ n B) - F(A)¥{B)\ : A G a(Y ),B G T n }. 

Let Y n = Y n - E(F ) and T n = ET=o Y t . 

^/E«=i X^ -2 "*^ 2 ^ < °°> i^-en T n /y/n converges in distribution to a Gaus- 
sian random vector = (T^t , • • • ,Too ) wtt -zero mean such that 

oo 

(3.3) E(T»rW) = E(F W F^' ) ) + ^ E(F (i) F^ + F^F«), 

n=l 

where the last series converges absolutely. 

The lemma is a combination of a one-dimensional CLT for mixing se- 
quences (cf. [10], page 425) with the multidimensional Cramer- Wold device 
(cf. [10], page 170). 

The next lemma is similar in spirit to [17], Lemma 2.1.10 and the proof 
of the latter works nearly verbatim. 

Lemma 3.2. Let the conditions of Theorem 1.3 hold. Then 

lim sup 1/n log Xn < 0, 

n— >oo 

where x n are defined by (3.2). 

Proof. Let 2) n denote the sequence (Xi)i>n, define To = inf{i G N : Uq^ = 
0}, and recall the definition of the mixing coefficients a(n) from (1.6). On 
one hand, using inequality (1.8), we obtain for any n G N, A £ o-(Yq) and 
B £ a(Yi-.i> n), 

F(Y G A, 2) n G B) > P(y G A r < n/2, 2)n G B) 

= Bp(P w (Yo € A^o < n/2)P w (2) n G B)) 

> P(y G A r < n/2) • P(2) n £ B) — 4o([n/2]) 

> P(y G A) • P(2) n £ B) — P(r > n/2) - 4a([n/2]), 
where [n/2] stands for the integer part of n/2. On the other hand, 

P(Y G A 2) n G B) < P(Y G A t < n/2, 2) n G B) + P(r > n/2) 

= Bp(B,(Y G A r < n/2)P w {% n G B)) +P(r > n/2) 
<P(Y G Ar <n/2)-P(2) n GB) 

+ 4a([n/2])+P(r >n/2) 
< P(Y G A) • P(2) n G B) + 4a([n/2]) + P(r > n/2). 
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Thus it remains to show that lim sup n ^ oc 1 /n log P(to > n) < 0. Using Cheby- 
shev's inequality in the third step and Jensen's inequality in the last one, 
we obtain 



P(r > n) = P > 1} < E F (C > 1) < E(||C/ ,n||l) 

\«=i / «=i 

= J Ep(||M n _ 1 M n _ 2 M io||i) 

< Ep(\\M n _i jK M n _2,K M , k /o,k||i), 

yielding the claim in virtue of assumption (Bl) since the matrix norm || • ||i 
and || • \\ K are equivalent. □ 

We next check that the second moment condition of Lemma 3.1 holds 
under the assumptions of Theorem 1.3. 

LEMMA 3.3. Suppose that assumption (Bl) of Theorem 1.3 is satisfied. 
ThenE(\\Y \\ K K )<oo. 



Proof. Recall the definition of the norm ||x||a;,K> % G 
i = 1, . . . , d we obtain, by using Minkowski's inequality, 

oo 

Epi " 



from (1.4). For 



E((Yi l) r)=Ep(\\Y^t) 



n=l 



<E P \ 



.n=l 



<( E[^(II<IM 1/K 



,n=l 



Conditioning on Uq n _i and using Minkowski's inequality again we obtain: 



\U< 



0,nl 



E E 

j=l m=l 



<E^tll^« 



r(i) 



K II O.Tl — 1 lllx),K' 



where the random vectors Ln, m are the same as in the branching equation 
(1.3). It follows by using induction that 



Ep(M 



<Ep(\\M n _ hK 



M),re/ ||re), 



and thus the claim of the lemma follows from assumption (Bl). □ 

Let T n = J27=o Yi and recall the vector p defined in (1.5). 

Lemma 3.4. Let Assumption 1.1 and conditions (Bl) and (B2) of The- 
orem 1.3 hold. Then the sequence n~ x l 2 (T n — n ■ p) converges in law to a 
Gaussian random vector T^ with zero mean. Moreover, if condition (B3*- 4 )) 

holds for some i G {1, . . . ,d} then E([T^] 2 ) > for that i. 
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Proof. The vectors Y n satisfy conditions of Lemma 3.1 in view of Lem- 
mas 3.2 and 3.3. Thus we only need to check that the ith component of the 
limit random vector is nondegenerate if condition (B3^) holds. 

Let Rn(w) = E UJ (Y n ). By (3.3) and using the fact that Yq and Y n are 
independent in a fixed environment, we obtain for any i = 1, . . . , d, 



E([T«] 2 )=E([Y C 



(0 



- p- n + 2 E E[(y (4) - P «)(Y« - P W)] 



:E([FJ 



n=l 
,(0l2 



WV 2 



(P W ) 



+ 2^i?p[( J Rj ) i) -p«)( J R«-p«)] 



n=l 



:=E([r «] a )-Bp([flf>] a ) + a*. 
By Lemma 3.1 the following series converges absolutely: 

oo oo 

£ E[(y w - 4*>)(yW - P W)] = £ ^[(4^ - pW)(i?« - p«)]. 

n=l n=l 

Moreover, it follows from Minkowski's inequality that 



E P ([R$ 



£p 



n 2\ 



X;(M n _i---M J ) (<) 

Ln=l 



< 



Ln=l 



and hence £'p([i?Q^] 2 ) < oo in view of (Bl). Therefore (see the last few lines 

in [5], page 198), a, = lim^ ^(Kpjflf " n " P^?) > °- Furthermore, 

if condition (B3«) holds, Y^ is a nondegenerate random variable under P u 
for P-a.e. environment uj. It follows that 



: j y -^ P ([/^Y)=^p[^(WT)-[^o w )n>o, 

completing the proof of the lemma. □ 



.(i)i2x 



^)m2i 



To complete the proof of Theorem 1.3 observe that 

n— 1 oo n— 1 oo 

^-^ = EE^< E E^M 



i=0 J=n 



i=— oo j=n 



and hence, in virtue of assumption (Bl), E(T n — S n ) is bounded by a vector in 
^ : E(Er=-oo ET=n Uij) = E(E,=-oo Ef o Uij) G M d . It follows that (T n - 
Sn)/y/n converges to in probability, yielding the claim of Theorem 1.3. 
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4. Recurrence behavior of the branching process. This section is de- 
voted to the proof of Theorem 1.5 and Theorem 1.6. The proof of the former 
is included in Section 4.1 while that of the latter is contained in Section 4.2. 



4.1. Proof of Theorem 1.5. First we make explicit the obvious link be- 
tween the structure of the limiting distribution and the recurrence properties 
of the sequence (Z n ) ne z + - For any environment cj, let 

(4.1) ttu(v) = P u (Z = v), veZ d + , 

and let tt(v) denote Ep(n w (v)) = P(Zo = v). According to Theorem 2.1, ir w 
is a probability distribution on Zl for P-almost every environment uj. In 
the following lemma we use the notation Z = (Z n ) n >Q an Z = (Z n ) n >o. 

Lemma 4.1. Let Assumption 1.1 hold. Then: 

(a) P(Z n = v i.o.) G {0, 1} for every v G lA. Moreover, ¥(Z n = v i.o.) = 1 
if and only if tt(v) > 0. 

(b) // tt(0) > then F(Z G -\Z = 0) = f n P^{Z G -)P(duj), where the 
probability measure P on the space of environments is defined by 

= ^7q^ ■ I n particular, the measures ¥(Z G •) and ¥(Z G -\Zq = 0) are 
equivalent if and only if P(7r w (0) > 0) = 1. 

Proof, (a) This is a direct consequence of Lemma 2.2. By part (c) 
of the lemma, P(Z n = v i.o.) = f(Z n = v i.o.), while part (b) implies that 
P(Z n = v i.o.) G {0,1} and furthermore, f(Z n = v i.o.) = 1 if and only if 
P(Zo = v) > (the latter follows from Poincare recurrence theorem in one 
direction and from the Borel-Cantelli lemma in the other one). 

(b) For any measurable set A C (Z^) z+ we have 

P(Z^nZ = 0)_ 1 :Ep[Pu{ z £Al z = 0)Pu{ z = 0)] 



P(Z = 0) 7T(0 

,T«(0) 



^(ZgA)- 



7T(0) 



□ 



It should be noted that P^^^O) > 0) G (0, 1) in a fairly common situation. 
Consider, for example, an environment which is a Markov chain in the state 
space {a, /3} and a process Z n such that ¥(Z\ = 0\coq = a) = while W(Z\ = 
0\luo = (3) = 1. However, using the fact that ^^(O) > i: L0 {v)P L0 {Z n = 0\Zq = 
v) for every v G the following simple criterion was in fact proved by Key 
in [14] (see the proof of the second part of [14], Theorem 3.3). 
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Lemma 4.2 ([14]). Let Assumption 1.1 hold and suppose in addition that 
for some uj G £1 there exists n G N suc/i i/iai P UJ {Z n = Q\Zq = 1) > 0, where 
1 = (1, . . . , 1) G K d . Then, TT$n w (0) > probability one. 

We note in passing that if the environment is an i.i.d. sequence and the 
law of the immigration q n is independent of the branching law p n -\ for all 
n G Z, then, integrating the equation irg^^v ) = X) u gz^ C^i =v\Zq = 

u), we obtain that 7r is a solution of the eigen- vector problem 7r = 7r^4 for 
the stochastic matrix A U)V = P(^i = i>|2o = u). 

Part (a) of Lemma 4.1 combined with Lemma 4.2 imply that under the 
conditions of the theorem, ¥(Z n = i.o.) = 1. The rest of the subsection is 
devoted to the proof that P(lim ri ,_ >00 u n /n = l/vr(0)) = 1 if (1.11) holds. 

Let P be the probability measure on the product space O x (Zl) z + de- 
fined by setting P(u G A, {Z n ) n&+ G B) = P(oj G A, (Z n ) neZ+ G B\Z = 0) 
for measurable sets ACQ and B C (Z^) z+ . The expectation with respect 

to P is denoted by E. For n G N let fi n = v n — v n -\ and consider the triples 
x n = {0 Un uj, fi n , Z\ Vn _ 1+ i Vn ]) , where the shift 6 is defined in (1.13) and the 
segments Zu^ in (1.12). The sequence (x ri ) n gpj is a Markov chain on the 

state space X = O x N x T d , where T d = Unez+ ( z i\{0}) n x {0}- We equip 
Yd with the discrete topology and denote by H the product cr-algebra in- 
duced on X. Transition kernel K of x n under P is given by 

K(u,i,x;A,j,y) =P w (Z [1J] =y)I e - jA (uj), 

where A G J-, i, j G N, x, y G T^, and wE(!. The initial law of x n is given by 

Q(A, i, x) := P(w EA,vi = i, Z [hui] = x) 

jFU^pL.i] =x)I -i A (o;)P(tL;), 

where the measure P is introduced in part (b) of the statement of Lemma 4.1. 

It follows from Lemmas 4.1 and 4.2 that if (1.11) holds, then P is equiv- 
alent to P and hence Q is equivalent to P. Moreover, applying Kac's recur- 
rence theorem ([10], page 348) to the stationary sequence Z n , we infer that 
E(/ii) = 1/tv(0). Therefore, the second part of Theorem 1.5 is implied by the 
following proposition. 

Proposition 4.3. Let Assumption 1.1 and condition (1.11) hold. Then 
the Markov chain (x n ) ng pj on (X, S) with initial distribution Q and transition 
kernel K is stationary and ergodic. 

Proof. The conclusion that the measure Q(-) = P(-|Zo = 0) is preserved 
under the action of the kernel K is standard (see, e.g., [15], Chapter II or 
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[16], Section 2.3) and follows from the fact that P is preserved under the shift 
on the sequence y n = (8 n u, Z n ). Specifically, for integers i,j such that i<j 
let Zuj] = (Zi, Zi+i, . .. ,Zj) and observe that for any i G N, x G Td, A G J 7 , 

Q(A, i, x) = ^nZo = 0, Z [1A = x, 9'u G A) 

= ^f(Z-i = °> ^[-i+1,0] = x, w G A). 
Therefore, using the notation ifQ(A, j, y) := f^K(ou;i,x.; A, j,y)Q(duj, i, x), 

KQ(Ai,y) 

dP(w) 



^2^2 p U z -i = °> z H+i,o] = x)P w (Z[u] = y)I fl -i j4 ( 



vr(O) 



= jf PUZ = 0, = y^^^M = Q(A, j,y). 

Thus Q is an invariant measure for the transition kernel K. 

Our next aim is to show that the sequence (x n ) n ^ is ergodic under P. The 
proof is a variation of the "environment viewed from the particle" argument 
adapted from the theory of random motion in random media (see, e.g., [4], 
Lecture 1 and references therein). Let A G S® N be an invariant set, that is, 
G" 1 ^ = A, P-a.s., where O is the usual shift on the space X N : (®X) n = x n+ \ 
for X = (x n ) n <=jq. It suffices to show that ¥(A) G {0, 1}. 

For x G X let ¥ x be the law of the Markov chain (x n ) n ^ with initial 
state x\ = x. Set h(x) = F X (A) and note that the sequence (/i(x n ))neN form 
a martingale in its canonical filtration. This follows from the following rep- 
resentation of h{x n ): 

E(Ia(A0|x O , Sl > ■■■■> X n)= E(I 6 -nA (X)\XQ, Xl, ...,X n ) 

(4-2) 

= h(x n ), P-a.s., 

where the first equality is due to the invariance property of the set A while 
the second one follows from the Markov property. The representation (4.2) 
yields, by Levy 0-1 law (cf. [10], page 263), that h(x n ) converges to Ia(X) 
as n approaches to infinity. This in turn implies that there exists a set fi£S 
such that h(x) = Ib(x), Q-a.s. (see the proof of (1.17) in [4]). 

It remains to show that Q(B) G {0,1}. By the martingale property of 
the sequence h(x n ), Ib(x) = KIb(x), Q-a.s. Since the transition kernel 
K(u,i,yi;-) does not depend of i and x, the function h should not depend 
on these two variables as well, that is, modulo a Q-null set, the set B has 
the form 

5 = {(w,i,x)GX:we(7} = CxNxT d 
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for some^C G T. 

Let K be the transition kernel of the Markov chain ui n = 9 Un u>, that is, 
for i G N and x G T^, 

K(u,A) = J2 Yl K(uj,i,x;A,j,y) = J2PUvi=j)h-JA(u)- 

Then 

jeNyeT d 

= Is(a;,i,x) = Ic(w), P-a.s., 

where the second equality follows from the fact that KIb = Is, Q-a.s. Thus, 
we have proved that 

Since P-a.s. the random variable Ic takes only two values either or 1, 
-FL(a*1 = 1) > and X)jeN-^>(Mi = i) = 1) we obtain that 

I c M = Ic7(^), P-a.s. 

That is, Ic(w) is an invariant function of the ergodic sequence (u; n ) ng z- From 
this it follows that P(C) G {0, 1} and hence Q(B) G {0, 1}. This completes 
the proof of the proposition. □ 

4.2. Proof of Theorem 1.6. First, we construct in an enlarged proba- 
bility space a representation of the law of Z n as a convex combination of 
an "external" multitype Galton-Watson process independent of the envi- 
ronment and a modified MBPIRE. The law P e of the external process is 
defined as a probability measure on (7~, Q) such that X n G lA. are i.i.d. with 

distribution e^ d+1 \ Ln,m G are i.i.d. independent of (X n ) nG z, and the 

law of Lnjm is e* for all n,m [the sequence {Z n ) n£ i + is defined as before by 
equation (1.3) with initial condition Zq = 0]. Roughly speaking, we attach 
to each particle a Bernoulli random variable independent of the environment 
in such a way that if the variable takes value 1 then the particle produces 
progeny according to the law P e . Thus, only if the random variable associ- 
ated with the particle takes value 0, the number of its progeny depends on 
the environment. 

We next show that for arbitrary large Iq G N there is a positive probability 
(that does not depend on i) that the branching process does not use the 
environment in time between vi and Ui + Iq and, moreover, i^+z = Vi + Iq. 
Due to the assumed mixing property of the environment, the branching 
process starts after occurrence of this event in the environment 9 Ui+l ouj which 
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is distributed "almost" according to P. The claim of Theorem 1.6 is derived 
then from this observation. The construction used in the proof is similar to 
that of L-safeguards for a RWRE introduced in [17], Section 3.1 and [7]. 

Let £ = (£„)„ e z+ and A = {Xn,m)nez+,meN,i=l,...,d be two collections of 
i.i.d. Bernoulli random variables on a probability space (E,£,Q e ) such that 
£ and A are independent of each other under Q £ and 

a(£n = l)=Q e (A« m = l)=£, 

Q 6 (£„ = 0) = Q(AW TO = 0) = l-e, 

where e is introduced in the conditions of the theorem. 

Further, let P denote the product measure on (fi x E,J- x £) whose 
marginal on (f2,.F) is P and that on (E,£) is Q E . The triple (w,£,A) serves 
as a new environment for the branching process (Z n ) ne z + defined as follows. 
We assume that for any realization (a>,£, A) of the environment, a law Pu>,£,\ 
on the underlying probability space (T,Q) is defined in such a way that: 

• Similarly to (1.3), (^ n )nGZ+ is a nonhomogeneous Markov chain that sat- 
isfies initial condition Zq = and branching equation 

d Z n +X n 

Zn+l = Yl E Z Sm forn>0, 

i=l m=l 

where 

• -Xn = I{ ?n= i}A" ni i + I{£ n=0 }X n , 2 , 



. fW _t r(») it rW 

i, n , m - A {A (*) m=1} ^n,m,l + 1 {A « m=0 } jL «,m,2> 

• each one of the four sequences (X n) i) neZ+ , (^„,2)nez + , 

(41,i)nez,mez + ,i=i,..,d and (^ i ) mi2 )nGZ,mGZ+,i=i,...,d consists of indepen- 
dent random variables, the sequences are mutually independent under 
Pu,,$,x, and 

(i) l \ (*)/ \ 

-'t^Av n,m,2 ~ V J ~ T" - "^ ' 

^,A(^n,l=^)=e^ +1) W, 

p /y x g^nW-gek » 

-rL),£,Al^n,2 = = , 

1 — £ 

where the constant e and the measures are defined in the statement of 
the theorem. 

Define the annealed probability measure 

(4.3) $ = P®P^,x = P®Qe®P^,\ 
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on the measurable space {£lxExT,J r ®£®Q). It is easy to see that the 
law of the process {Z n )nez+ under P coincides with that of {Z n ) n& i + under 
P, while its law under Q e (8) Pw£,\ is the same as that of (Z n ) ne z+ under 
P u . Therefore, part (a) of Theorem 1.6 is contained in the following lemma. 
Part (b) of the theorem follows from the lemma in virtue of Lemma 4.1(b) 
and Proposition 4.3. 

Lemma 4.4. Let the conditions of Theorem 1.6 hold. Then 

sup sup sup {P(A|B)-P(A)}^ n ^ oo 0, 

m ^Aeg""+™ Bee„ m ,P(B)>o 

where we denote Q n = a{Zi :i>n) and Q n = a{Zi,oJi :i <n). 

Proof. With a slight abuse of notation let 

uq = and u n = inf{i > v n -\ '■ %i = 0}. 

Fix any Iq G N. By the conditions of the theorem, P e (f/ = lo) > 0. This 
implies that for some integer R, 



(4.4) P e 1% = Z and max \\X t i |h < R > 0. 

V. o<i<; — i 11:11 ) 

Define the sequence of independent under P events (A)iez+ by setting 
Hi = U Vi+t = 1; >^-+t,m = 1 for * = °) • • • > h - M = 1, • • • , d, 

m = l,...,R;u i+ i =Ui + l ; max < R }■ 

Vi<3<v i+ i Q -l ) 

Note that on the event Hi, the branching process does not use the envi- 
ronment during Iq = fj+/ — Vi steps after z/j. It follows from (4.4) that the 
probability of Hi, P(A), is positive and does not depend on i. 

For any n G N define m n (lo) = min <j< n _/ {i : Hi occurs}, with the con- 
vention that the minimum over an empty set is infinity. For any i <n — Iq 
and t £ N, let A,t = {m n (lo) = i,Vi = t} and 

A = |J A,t = {m n (l) = i} 
ten 

and 

n— lo 

D= (J D t = {m n (l ) <oo}. 

i=0 

The events A are disjoint and can be represented in the form 

A = # c n---n#tin A- 
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Since F(Hi) =:p>0, we obtain that F(D) = EIU' "^ 1 ~ pYp ^n-*oo 0. 

Let be the shift by v\ on the sequence Z = (Z n ) ne z+, that is, (@ n Z)i = 
Z Vn+ i. For t G Z denote P* = P(g> Q e <8> ■Pe* W) e Aj where 6* is the shift operator 
defined in (1.13). Using inequality (1.14), we obtain that for any events 

n—lo 

F(A\B) <F(A, D\B) + F{D C \B) = J2 J2 D i,t\ B ) + 

i=0 t 

n— io 

= £ £p(A|A l t,5)P(A 1 t) + P(£ c ) 

i=0 i 

< EE( ?+lD ( e " ! " i0A ) + M*o))P(A,t) + P(£ C ) 
i=0 t 

n—lo 

< P(e~ i_i M)P(A) + 2vp(Z ) + p( j d c ). 

i=0 

Similarly, 

ra— Jo 

P(A|B)>P(A, J D|S)=5^53P(A,A,*|-B) 

i=0 t 

1=0 * 
n— Zo 

> £ ^ ( pt+Zo (e -i-Zo A ) _ 2p(Z ))P(A,0 
i=0 a; 

n— io 

> £ p(e- i - i M)P(A) - 2^(z ). 

i=0 

It follows that |P(il|J3) - F(A)\ < 4((f(l ) + F(D C )). That is, for n large 
enough we have |P(A|S) — P(-A)| < 8ip(lo), proving the lemma since Zo is an 
arbitrary integer. □ 
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